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Abstract 

We present a theoretical study of wetting phenomena and interactions between liquid-vapor 
interfaces based on the local density functional theory. The focus is mostly on the impact of long- 
range van der Waals interactions both within the fluid and between the fluid and substrate. For 
the latter, we consider two models - hard wall and soft wall approximations - differing by the role 
of steric effects and leading to a qualitatively different character of phase transitions. We compute 
numerically the disjoining and conjoining potentials (which are important dynamically for spreading, 
spinodal dewetting, and coarsening in thin films, as well as resolution of interfacial singularities), 
and loci of intermediate and complete wetting transitions as functions of the Hamaker constant and 
temperature. We find that taking into account short-range interactions is essential for the description 
of wetting transitions in the soft wall limit. We derive an analytical form of the disjoining potential 
and analyze it in the context of the complete, frustrated and partial wetting. 

1 Introduction 

Recent studies of dynamic behavior of thin fluid films and three-phase contact lines showed the crucial 
role of precursor layers in contact line dynamics and droplet interaction in the course of spinodal dewet- 
ting and coarsening ^ |2 CU ■ While the role of precursor layers in wetting and spreading phenomena 
had been discussed for a long time 0], only recently it was realized that the motion of contact lines is 
strongly influenced by the precursor layers even when the fluid is not completely wetting ^ [S] . 

Intermediate wetting transitions leading to the formation of mesoscopic fluid layers, as distinct from 
complete wetting when a macroscopic layer is formed, were discovered long time ago jSJ [7J |51 |5] , and 
observed recently experimentally on solid ^U] an d liquid ^2 021 El [HI El substrates. Both complete 
and intermediate wetting transitions are assumed to incorporate both short and long-range interactions. 
However, a well formulated theoretical study which justifies the individual impact of such interactions 
is still missing, though it is recognized as a prime theoretical subject of interest in view of recent 
experimental observations and their partial agreement with available phenomenological models 16 . 

The mesoscopic dynamic diffuse interface theory ^ [z] uses the wide separation between the width 
of mesoscopic layers and the scale of hydrodynamic motion to reduce the hydrodynamic equations to 
a generalized Cahn-Hilliard equation where the disjoining potential of the fluid layer |17j serves as 
the driving force of the fluid motion along the film. A similar approach applied to a mesoscopic layer 
separating two identical bulk fluid phases can be used for the description of coalescence dynamics, and 
has been recently discussed in connection to resolving an interfacial singularity in viscous flow |18j . 

The cited works used model expressions for the disjoining potential obtained in the sharp interface 
approximation, while emphasizing the essential role of the diffuse character of the vapor-liquid inter- 
face (known already to van der Waals JH1) f° r resolving the interfacial singularities. The dynamic 
importance of the disjoining potential in mesoscopic layers necessitates, however, precise computation 
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linking it to molecular interactions in the fluid and between the fluid and substrate. The adequate 
tool for these computations is density functional theory (DFT) [201; further reduction to a local (van 
der Waals-Landau-Cahn) theory through small gradient expansion is invalid when the interaction is 
long-range (i.e. has a power decay Numerical DFT computations were carried out in the con- 

text of phenomenological phase transition studies [51l22| without emphasis on computation of disjoining 
potential. 

In this communication, we explore this problem anew using the simplest long-range interaction po- 
tential, which allows for a direct link between the wetting properties and thermodynamics of a van der 
Waals fluid. We pay special attention to the character of interactions in the vicinity of the substrate, 
distinguishing between "hard wall" and "soft wall" limits. The first version used in earlier computa- 
tions [22| exaggerates steric density depletion near the substrate, and the opposite limit may be more 
relevant for rough solid interfaces, as well as for films on liquid substrates. The character of wetting 
transitions interactions, as well as properties of the precursor layer, are very much different in both 
limits. Realization of the latter reflects the droplet roll-up behavior under different substrate surfaces; 
for which the contact angle passes tt, as was observed in recent experiments [23 E]. It is found that the 
analysis in the soft wall approximation based solely on van der Waals disagrees with results of recent 
experiments on binary fluid systems ^2 1121 1131 1141 IT5| . We included therefore a weak dependence on 
short-range interactions; the expression for the disjoining potential modified in this way was found to 
be in a good qualitative agreement with the experiment. 



2 Density functional equations 

Our starting point is the free energy functional written in the local density functional approximation |25| 
as 

T= j p(r)/[p(r)]d 3 r+i J p(r)d 3 r^(7(r')[p(r + r')-p(r)]dV, (1) 

where f(p) is free energy per particle of a homogeneous fluid and U(r') is an isotropic pair interaction 
kernel with a short-scale cut-off d. The functional Q is written in the form |26j separating the contri- 
bution of density inhomogeneities, expressed by the last term vanishing in a homogeneous fluid, but is 
equivalent to an alternative form used in earlier works |21| . 

The chemical potential p = 8J- /Sp enters the respective Euler-Lagrange equation obtained by min- 
imizing the grand ensemble thermodynamic potential $ = T — p J pd 3 r, which defines the equilibrium 
density distribution p(r): 

g(p)-p+ ( U{r')[p(v + r') - p(r)] d V = 0, (2) 



where g(p) = d[pf(p)]/dp. The function F(p) = p[.f(p) — p] should have two minima p^ corresponding 
to two stable uniform equilibrium states of higher and lower density (liquid and vapor). 

A simple example of long-range potential is the modified Lennard- Jones potential with hard-core 
repulsion: 

U={ ~ Cwr ^ & ' r> i , (3) 
|_ oo at r < d w 

where d is the nominal hard-core molecular diameter. The interaction kernel U(r) gives the free energy 
density of a homogeneous van der Waals fluid p)| 

f(p,T) = T In -JL^-ap, (4) 

where T is temperature, b = |7T<i 3 is the excluded volume and 

2 n 2-kCw 



a = -2ir I U(rydr=-^-. (5) 



Equilibrium between the two homogeneous states, p — p is fixed by the Maxwell condition 

Pof(pp) - Pof(pp) /o\ 

Mo — r z , (oj 

Po ~Po 
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which defines, together with po = g(p ), the equilibrium chemical potential p = po and both equilibrium 
densities. 

The equation for the density distribution near a flat boundary normal to the z axis is obtained by 
assuming p to be constant in each lateral plane and integrating Eq. Q in the lateral directions. This 
yields the free energy per unit area, or surface tension 

oo oo 

p{z)[f{p)-p]dz+ l - J p{z)Az J Q(0[p(z + 0-p(z)}dC. (7) 

— oo — oo 

The interfacial energy is contributed both by deviations from the equilibrium density levels in the 
transitional region and by the distortion energy localized there. The ID interaction kernel Q(z) lumps 
intermolecular interaction between the layers z = const. It is computed by lateral integration using as 
an integration variable the squared distance q = r 2 = £ 2 + z 2 , where £ is radial distance in the lateral 
plane. Taking note that the lower integration limit for q is qo = z 2 at \z\ > d, q = d 2 at \z\ < d, we 
compute 

„oo ( —\-kCwz~^ at \z\ > d 

Q(z) = -7rC w q- 3 dq=\ ' (8) 
J io [ -±irC w d- 4 at \z\ < d. 

The respective ID Euler-Lagrange equation, replacing Eq. (J2J), is 

Q(0[p(z + Q -p(z)] d( = 0. (9) 

-oo 

This equation can be rewritten in a dimcnsionless form 

n poo 

g(p) -H+jP QKMz + C) - p(z)\ dC = 0, (10) 

J — oo 

where 

g(p) = jL- - In (~ - l\- 2 pp. (11) 

Here the length is scaled by the nominal molecular diameter d, the density by £> -1 , and the chemical 
potential by T; the interaction kernel is Q(z) — —z~ 4 at \z\ > 1, Q(z) = —1 at \z\ < 1, and the only 
remaining dimensionless parameter is the rescaled inverse temperature [3 = a/(bT). 

An example of a density profile obtained by solving numerically Eq. I|10(l is shown in Fig. ^ The 
density tail asymptotics can be estimated by considering a location far removed from the interface 
placed at the origin [\z\ ^> 1 in the dimensionless units of Eq. I|10|l ] where a sharp interface limit can be 
implemented. The density is presented as p = p^ + p, where p/p ~ l/|z| 3 *C 1 . Inserting this in (| 1 1 1) 
and linearizing around p = p^ , we see that the densities inside the integral are well approximated in the 
leading order by the two limiting constants, which is equivalent to the sharp interface limit. For example, 
for the vapor tail at z > 0, \z\ > 1 we have p(z) = p$ and p(z + () — p^ for Q > \z\, p(z + C) = Po f° r 
C < \z\. Thus, we obtain 



P(Pa - Pa 



This is in good agreement to the numerical solution, as seen in the inset of Fig. One can check 
a posteriori using this expression that the contribution to the integral of neighboring locations with 
|£| = O(l) is of a higher order oc |z|~ 5 and therefore can be neglected. 



3 Interacting fluid-vapor interfaces 

If two flat fluid-vapor interfaces are in close proximity, the equilibrium chemical potential and surface 
tension become dependent on their mutual separation h. This phenomenon is relevant for processes of 
topology change, e.g. droplet coalescence. The corrections due to proximity of interfaces in the case when 
a gas layer separates two identical semi-infinite bulk liquids can be obtained by finding a homoclinic 
solution of Eq. (|l()(l satisfying the asymptotic conditions p(±oo) = p + . A stationary solution of such 
kind exists at certain values of p shifted downwards from the Maxwell construction level p,Q. The shift 
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Figure 1: The density profile of the liquid- vapor interface at p = 9 obtained by numerical solution 
of IjlUI) . The inset shows the vapor-phase tail of the numerical solution (solid line) compared to the 
asymptotic form l|12|l depicted by the dashed line. 



A/i c = /io — /i corresponds to the conjoining potential expressing the interaction of two identical flat 
interfaces. 

A rough but qualitatively correct approximation can be obtained by computing molecular interactions 
between two sharp interfaces |18| . The shift of chemical potential necessary to keep two interfaces 
separated by a distance h in the state of equilibrium is determined in this approximation by the decrement 
of the integral in Eq. I|1(J|I due to replacing gas by liquid at z > h. For the purpose of a qualitative 
estimate, the sharp-interface computation valid at h 3> 1 can be extended also to small separations, 
(h < 1) to find QH 

{§ - h at h<l, 
(13) 
i/i~ 3 at h>l, 

Equilibrium of a layer between two interfaces is unstable; the instability is, however, very weak when 
separation of interfaces is large compared to the molecular scale. Localized small perturbations decay 
under these conditions due to surface tension, and a large disturbance is needed to create a critical 
nucleus initiating the topology change. 

A precise dependence is obtained by solving numerically Eq. ()10(l . The solution is found by fixing 
some trial value of p and solving Eq. (|10J) iteratively to find a stationary profile p(z) at this value. The 
nominal gap width is defined as 




The computation results for (3 — 9 are shown by dots in Fig. [5] The curve A/i c (/i) expressing this 
dependence well fits the computational results shifted by a certain value h+, equal to « 1.39 in this 
particular computation. A shift is necessary because, while the separation in Eq. (|13(l can be formally 
arbitrarily small, no stationary solution of Eq. IjlUf) can exist below a certain value of h which corresponds 
to a critical size required for nucleation of a critical ID "bubble" . The applied shift equals to the width 
of this "bubble" computed according to Eq. I|14|) . 

4 Fluid-substrate interactions 
4.1 Disjoining potential 

In the proximity of a substrate surface, the additional term in the free energy integral £Q| is 

T.= l P (v) d 3 r I U s (\r - r'\)p s (r') dV , (15) 
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Figure 2: The dependence of the dimensionless conjoining potential Ap c on separation h. The dots 
denote the results of ID density functional computation with a shift of h* = 1.39 adjusted to fit Eq. I|13fl . 
as shown by the solid line. Inset: the tail of the numerical dependence Ap c (h) compared to the h~ 3 
decay (dashed line). 



where U s is the attractive part of the fluid-substrate interaction potential, p s is the substrate density, 
and J means that the integration is carried over the volume occupied by the substrate; all other integrals 
here and in are now restricted to the volume occupied by the fluid. 

In the following, we shall consider a flat interface parallel to the substrate surface z — 0, and 
suppose that liquid-substrate interactions are also of the van der Waals type with a modified constant 
Cs = a s Cw- Then the free energy per unit area is expressed, after some rearrangements, as 



p(z)\f(p) + Mz) 



dz + 1 / ,,[;),!: / ()(: - 0[/'k >-/;(:>] <k. (1<>) 



The first term contains the same local part as in Eq. Q complemented by the liquid-substrate interaction 
energy. The latter is computed by integrating the attracting part of the fluid-fluid and fluid-substrate 
interaction energy laterally as in Eq. JHJl and represents the shift of energy compared to the unbounded 
fluid. The term p(z)/2 compensates lost fluid-fluid interactions in the substrate domain which are 
included in the homogeneous part f(p). 

Computation of the function ipi( z ) depends on steric constraints imposed upon fluid molecules in the 
vicinity of the substrate. We consider two limiting cases, as demonstrated in Fig. 

(i) Soft wall approximation (SWA) which allows for fluid molecules to penetrate up to the surface of 

the substrate. 

(ii) Hard wall approximation (HWA) which imposes steric constraints preventing the centers of fluid 

molecules from approaching the substrate at distances shorter than the hard-core cutoff d. 

To be definite, we place the origin of the coordinate system on the centerline of the first row of the 
substrate atoms. In the soft wall approximation, the computation yields 

/0 poo />0 

dC / q- 3 dq= / Q(z-Qd(, (17) 
-oo J qo J— oo 

where the integration limit is go = ( z — C) 2 at \z — CI > d, q = d 2 at \z — £| < d. The result is 

!-\nC w z~ 3 at \z\ > d 

(18) 
-nCwd- 3 (| - at \z\ < d. 

The dimensionless Euler-Lagrange equation derived from Eq. (|16[) reads 

oo 

g(p) - p + \f3Uz) [ P + (X + 1) - P(z)] + |8 / Q(C ' z)[p(0 - P {z)\ dC = (19) 
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Figure 3: Schematic cartoon representation of soft (left) and hard (right) wall approximations, 
circles represent hard sphere fluid molecules and the horizontal coordinate, in this setup is z. 
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Figure 4: The function if)i(z) for soft (solid line) and hard (dashed line) wall approximations. 



where \ = &sPs/ P + — 1 is the dimensionless Hamaker constant, and ipi( z ) = — z 3 /3 at z > 1, ipi( z ) = 
z - 4/3 at z < 1. 

The hard wall approximation introduces, in effect, a void layer with the thickness equal to the 
hard-core cutoff, hence i/i[(z) = at < z < 1 [22] • If the closest allowed position of the centers 
of fluid molecules, is taken as the origin, then we define a hard wall fluid-substrate interaction as 

as seen in Fig [J] This shift significantly changes the equilibrium solution and 
the character of the wetting transition. 

The equilibrium chemical potential is shifted from the Maxwell construction, /i = /io in the proximity 
of the substrate surface. The shift A fid = fj, — [j,q, called disjoining potential fJ7j, can be defined as 



1 d ls 
Po - Po dh 



(20) 



where h is the nominal distance between gas-liquid and liquid-substrate interfaces. The latter is defined, 
analogous to Eq. I)14p. as 

1 f°° 

h / (p-p-)dz. (21) 



P Jo 



4.2 Equilibrium solutions 

The sign of x distinguishes a perfectly wetting fluid at x > and a "nonwetting" fluid forming a 
finite contact angle at x < 0. In the latter case, there are two branches of equilibrium solutions of 
Eq. (|19|l parametrized by the disjoining potential A/i^. The stable branch with small h is characterized 
by a monotonic density profile and corresponds to the vapor phase thickening to a relatively dense 
adsorbate or precursor layer near the substrate. The unstable branch with larger h is characterized 
by a non-monotonic density profile and corresponds to a liquid layer with a slightly depleted density 
near the substrate. Instability is characteristic to any layer of a nonwetting fluid, but it is very weak 
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Figure 5: Equilibrium solutions of Eq. (|19fl in the soft wall approximation (SWA), (a) The depen- 
dence of the dimcnsionless disjoining potential Afid on separation h. The solid/dashed lines denote 
stable/unstable solutions, respectively. The dotted line in the inset depicts the h~ 3 decay, (b) Coexist- 
ing density profiles at Afid = 2.2 • 1CP 3 . Stable (solid line) and unstable (dashed line) profiles correspond 
to S and U in (a) respectively. Parameters: \ — —0.05, (3 = 7. 



when separation of the vapor-liquid and liquid-substrate interfaces is large compared to the molecular 
scale. The contact angle can be expressed through the disjoining potential for the case |x| <C 1 when a 
small-angle approximation is valid |26| : 



0=\—r^d&h, (22) 
V 7 Jh 

where h = ho is the precursor layer thickness defined by the condition /u(/io) = fio or Afj,d(ho) = 0. 

Figures EI a) and EJa) present typical equilibrium curves Afj,d(h), respectively, for the SWA and 
HWA. Examples of the corresponding density profiles are shown in Fig.JSJb) and[|Jb). One should keep 
in mind that in the hard wall approximation, h should be rescaled back so that the density profiles 
start at ft. = 1. Unlike in the soft wall approximation, in the hard wall case, the void between the 
substrate and the liquid density encourages density depletions [see inset of Fig. E£b)]. All solutions 
exhibit a A/i^ ~ h~ 3 tail at large h, in agreement to calculations performed in the sharp interface 
approximation |26|. Oscillatory density tails cannot appear in our model, unlike more sophisticated 
nonlocal DFT computations |271 122|. Although the curves A/j,d(h) are qualitatively similar to those 
obtained in the sharp interface approximation |26| . the quantitative distinctions strongly influence the 
character of the wetting transition, as will be emphasized in the following. 

4.3 Comparison between HWA and SWA 

Investigating numerically the two above substrate-liquid interaction models we can distinguish between 
two main differences in the interaction properties: the emergence of "microscopic" solutions, identified 
with nanoscale precursor layers and transition to layers of mesoscopic or macroscopic thickness. The 
three classes of films correspond, respectively to h ~ 0(1), h ~ O(10), h > O(10 2 ) (measured on the 
molecular scale d). 

4.3.1 Precursor layer 

Stable equilibrium solutions with a finite thickness h — ho, which correspond to a microscopic precursor 
layer, exist at the liquid- vapor equilibrium chemical potential /i = /io- An example of the dependence of 
the dimensionless precursor layer thickness on the inverse temperature /3 is shown in Fig. Ufa). One can 
see a strong difference between SWA and HWA results. In the HWA computation, the value of ho at 
the chosen value of the Hamaker constant is much less than unity, unless near the critical temperature 
P = P c ~ 3.37, so that one can speak of a dilute adsorption layer rather than of a proper precursor. 
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Figure 6: Equilibrium solutions of Eq. (|19fl in the hard wall approximation (HWA). (a) The depen- 
dence of the dimcnsionless disjoining potential A/i^ on separation h. The solid/dashed lines denote 
stable/unstable solutions, respectively. The dotted line in the inset depicts the h~ 3 decay, (b) Coexist- 
ing density profiles at A/i^ = 7.7 • 10~ 3 . Stable (solid line) and unstable (dashed line) profiles correspond 
to S and U in (a) respectively. Parameters: x = —0.05, (5 = 1. 



This difference stems from an effective increase of the absolute value of the Hamaker constant due to 
the presence of a dilute layer of steric origin present in HWA. 

On the other hand, by fixing the value of (3 and varying x, one finds that a precursor layer may exist 
only above a critical value x > Xc i as shown in Fig.[?Ib). When the layer thickness is defined by the 
integral expression Q21[l. this transition loses a qualitative character, and the value xZ can be defined as a 
point where Afid(ho) = [see Fig.lHJa)]. This happens at x = XciP) = — (1 — P~ I ' P + )> so that the lower 
limit is identical for both models (see Fig. |5J . According to the integral formula, negative values are 
possible, and may appear when fluid-substrate interactions are so weak that the fluid is nonwetting even 
at vapor densities. Moreover, at x slightly above Xc ■> the topology of the curves Afid(h) for the HWA 
model changes: the curve becomes discontinuous, and the microscopic and macroscopic branches of the 
curve separate, as seen in Fig. [Hta) • The discontinuity is explained by the absence of the vapor- fluid 
coexistence above some critical value of the chemical potential as shown in the inset of Fig. Ufa) . The 
sequence of density profiles in the vicinity of x = Xc 1S shown in Fig. |8|V). 

4.3.2 Wetting transitions 

A precursor film exists within the interval xZ < X < xt ■ The latter limit corresponds to the wetting 
transition, as presented in Fig. |7fb). One can see a qualitative difference between the SWA and the 
HWA models. For SWA, the precursor thickness grows continuously as x increases, and diverges at 
X = xt = 0, indicating a second-order transition to complete wetting. For HWA, the precursor layer 
approaches a finite thickness at finite x = Xc > indicating a first order transition to wetting. The locus 
of the wetting transition x = xtifi) is determined by numerical computation (see Fig. EI> • 

The two limiting cases of SWA and HWA correspond to compliant (fluid or rough solid) or molecularly 
smooth solid substrates, respectively. Recently the transitions of both kind (first and second order) were 
found in experiments on liquid substrates |151 . The quantitative disagreement between the SWA 
results discussed here and in the later experiments is due to the absence of short range interactions in 
our computations. A more complete picture of wetting transitions in SWA arises when, in addition to 
long-range interactions, short-range forces between the fluid molecules and the substrate are taken into 
account ^3 EH- While concentrating in the following on short-range interactions in SWA, we refer the 
reader to Ref. [22] which describes the impact of short-range interactions in HWA. 

4.4 Short-range interactions in SWA 

Computation of repulsive short-range interactions in the soft wall approximation takes into account polar 
interactions. In the simplest description, the interaction kernel for short-range forces has an exponential 
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Figure 7: The dependence of the dimensionless precursor layer thickness on (a) the inverse temperature 
(3 for a fixed value of the dimensionless Hamaker constant x = —0.05; and (b) the Hamaker constant x 
at (3 = 4. The solid/dashed lines represent the soft/hard wall approximations, respectively. In (b) Xc 
corresponds to the emergence threshold of the precursor layer, identical in both SWA and HWA models. 
The precursor thickness at the wetting transition threshold x — xt i s finite in HWA, indicating a first 
order transition, while in SWA ho — > oo at x ~ > xt = 0, indicating a second order transition. 




Figure 8: Equilibrium solutions of Eq. (|19|l in the hard wall approximation (HWA) in the vicinity 
X = Xc = —0.833. (a) The dependence of the dimensionless disjoining potential A/^ on separation h at 
j3 = 4 and different values of the Hamaker constant, from right to left, x = —0.6, —0.8, Xc ■ The dashed 
line represents the upper limit of the vapor-liquid coexistence. The coexistence range is demonstrated 
in the inset. The dashed line denotes the critical shift of chemical potential at which the curve A/i^/i) 
becomes discontinuous, (b) Density profiles corresponding to the arrows in (a). 
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Figure 9: The dependence of the critical values of the Hamaker constant Xc on inverse temperature f3 
at (J, = fio- The left branch [xc(P) < 0] is the limit \I of the emergence of the precursor layer identical 
for both SWA and HWA models. The right branch [xc({3) > 0] denotes the wetting transition to for 
HWA. The wetting transition for SWA occurs at \ c — independently of f3. 



o 



i> s = 7tC p / dC / e~ X ^ dq , (23) 



decay QSII3Q: 



where A > 0. Some other forms can be also used here [3T]. Unlike the computation of ipi, no cutoff is 
required, and the lower integration limit is go = i z — Q 2 ■ This yields 

^ = HZ^k( 2 + Az) e - A *. (24) 
The modified dimensionless Euler-Lagrange equation reads 

oo 

g{p) ^s(z) + \Wi{z) [p+( X + 1) - p(z)] Q(C - zMO ~ P(*)\ dC = (25) 

o 

where for simplicity we set A = d^ 1 , i> s {z) = (2 + z) exp(— z) at z > and 77 = C p b 2 /a ~ C p /Cw- An 
analytical form of the disjoining potential can be obtained exploiting the sharp interface approximation 
for both rj and % <C !• Adapting the method discussed in 0121] we find 



A P 



UW2 + h)e-»-p+& P+(] 1 



/l 3 8 5 '(/7+)fc 6 



(26) 



The derived analytical form of the disjoining potential (|26l) is presented in Fig. I1UI A comparison with 
direct numerical integration of 125|) shows a good agreement already at h < O(10) values. 

A qualitative nature of wetting transitions can be captured via a simple analysis of the disjoining 
potential l|26ll at the equilibrium chemical potential fi = /iq and thus A/i c i = 0. We distinguish between 
three characteristic types of behavior which are attributed to the three corresponding wetting regions: 

Complete wetting - there is a single diverging solution (h — > 00) at % > 0, A/i^ = 0, which corre- 
sponds to a macroscopic layer. 

Frustrated complete wetting - three solutions exist for x > 0, two of which are finite and the third 
one is diverging. The two finite solutions have distinct scales, molecular ("microscopic") and 
mesoscopic. 

Partial wetting - two solutions exist for x < 0, where one is finite on a molecular scale and the other 
one diverges. 
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Figure 10: The dependence of the dimensionless disjoining potential A/^ on separation h. The solid 
line denotes the sharp interface approximation form [Eq. (|26|l ] while the direct integration of Ij25(l is 
demonstrated by the dashed line. Parameters: = 9, r\ = 0.01 and \ = 0.3. 



The separating boundary between the partial and frustrated complete wetting depends solely on the 
sign change of the Hamaker constant, x- The transition from frustrated to complete wetting occurs at 
a certain critical value of n = r\ c when the maximum of the curve A^id(h) crosses zero; this happens at 
a critical thickness h = h c defined as 

<0. (27) 

h=h c 

The dependence A[id(h) at r/ = r\ c is shown by the upper curve in Fig. 1111 This corresponds to a 
discontinuous (first order) intermediate wetting transition from microscopic to mesocopic films. The 
lower curve with an inflection point in Fig. II II corresponds to a continuous (second order) intermediate 
wetting transition, which is observed at a certain shifted value of chemical potential A/i^ < 0. 

Figure IT2l summarizes these three transitions in the x ~ V plane. The critical point x = V = marks 
the critical end point |28| and corresponds to the sign change of the Hamaker constant; all three above 
wetting regions converge at the critical end point. Changing the inverse temperature (3 (for constant 
values of rf) does not change this picture in a qualitative way. 
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Figure 11: The dependence of the dimensionless precursor layer thickness on the disjoining potential 
according to Eq. |2E1 The lower line (rj = 0.00344) corresponds to the continuous intermediate wetting 
transition. The upper line (rj = r\ c = 0.00407) corresponds to onset of the frustrated complete wetting. 
The solid/dashed lines represent stable/unstable solutions, respectively. The Parameters: (3 — 9 and 
X = 0.4. 
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Figure 12: The boundaries of wetting regions (solid line) in the parametric plane x ~ V f° r = 9 and 
at the equilibrium chemical potential A/id = 0. The dashed line depicts the continuous intermediate 
wetting transition at A/i<j < 0. 



5 Discussion 

In the current study we have presented an analysis of wetting phenomena and interactions between liquid- 
vapor interfaces at microscopic, mesoscopic and macroscopic distances based on a relatively simple local 
density functional theory (DFT). The study is centered on computation of disjoining and conjoining 
potentials (which are important for dynamic computations). Focusing on the role of long-range van der 
Walls interactions, we considered two models - hard wall and soft wall approximations - differing by the 
role of steric effects. The first difference corresponds to a precursor layer thickness and the quantitative 
form of the disjoining potential. Under similar external conditions we found that droplets on smooth 
surfaces will roll-up before the ones on complaint surfaces |^] due to the steric forces. The distinct 
wall models also lead to a qualitatively different character of wetting transitions. A first order transition 
from partial to complete wetting occurs in the hard wall model only, turning to a second order transition 
in the soft wall model. 

To capture a more realistic description of wetting transitions, we included a weak dependence on 
short-range polar interactions. Analytical derivation of the disjoining potential was carried out exploiting 
the sharp interface approximation, which showed a good agreement with a direct numerical solution of 
DFT Euler-Lagrange integral equations. We find that the qualitative nature of the sequence of wetting 
transitions stems from a competition between short- and long-scale interactions, best seen in the (x, rj) 
parametric plane, while all other parameters, like temperature or short-range decay range A have 
only a quantitative effect. The role of repulsive short-range interactions is analogous to steric effects, 
which are emphasized in the hard wall approximation. Thus, the presented simple model reproduces 
all major types of wetting transitions (between partial to complete/frustrating wetting and thin/thick 
precursor in coexistence with bulk fluid) of |14l 1151 ITfil I28j as summarized in Fie. 1121 
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